NACA TN 3815 09fO~ 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



TECHNICAL NOTE 3815 


ON SLENDER-BODY THEORY AND THE AREA RULE 

AT TRANSONIC SPEEDS 

By Keith C. Harder and E. B. KLunker 

Langley Aeronautical Laboratory 
Langley Field, Va. 




Washington 
November 1956 


AFMDC 


TECVlI'T't 


OBSASiY 


- — J M 

\i 1 x 

iyu 


tech library kafb, nm 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTI 


T 


TECHNICAL NOTE J8l5 

ON SLENDER-BODY THEORY AND THE AREA RULE 
AT TRANSONIC SPEEDS 1 


TECH LIBRARY KAFB, NM 



mu ni:i mu Jin hi 

□ 0fefa70? 


By Keith C. Harder and E. B. Klunker 


SUMMARY 


The basic ideas of the slender-body approximation have been applied 
to the nonlinear transonic -flow equation for the velocity potential in 
order to obtain some of the essential features of slender-body theory 
at transonic speeds . The results of the investigation are presented 
from a unified point of view which demonstrates the similarity of slender- 
body solutions in the various Mach number ranges. The primary difference 
between the results in the different flow regimes is represented by a cer- 
tain function which is dependent upon the body area distribution and the 
stream Mach number. The transonic area rule and some conditions con- 
cerning its validity follow from the analysis. 


INTRODUCTION 


Slender-body theory originated with Monk's paper (ref. 1) in 1924 
in which the forces on slender airships were calculated for low-speed 
flight. In 1938 Tsien (ref. 2) pointed out that Hunk's airship theory 
also applied to the flow past inclined pointed bodies at supersonic 
speeds. The subject gained new importance in 1946 with the appearance 
of Jones's paper (ref. 3 ) in which it was shown that the basic ideas of 
the slender-body approximation could be used to calculate the forces on 
slender lifting wings at both subsonic and supersonic speeds provided 
that proper account was taken of trailing -vortex sheets . Since Jones ' s 
paper, the subject has received wide treatment. In an important paper 
in 1949# Ward (ref. 4) developed a general unifying theory for the flow 
past smooth slender pointed bodies at supersonic speeds . This theory 
contains as special cases the lifting planar wings of Jones and the slen- 
der no nli fting bodies treated by Von Karman (ref. 5)* The corresponding 
problem at subsonic speeds has been e xamin ed by Adams and Sears (ref. 6) 
who also extended the slender-body concepts to shapes which are "not so 
slender." Lighthill (ref. 7) has given a method for calculating the 


Supersedes recently declassified NACA Research Memorandum L54A29a 
by Keith C. Harder and E. B. Klunker, 1954. 
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flow past bodies with discontinuities in slope. Keune (ref. 8) has 
developed solutions for slender wings with thickness, and various lifting 
configurations have been treated by Heaslet, Spreiter, Lomax, Ribner, and 
others (refs. 9 to 13)* 

The Blender-body theory presented in references 2 to 13 has been 
based upon the linearized equation for the velocity potential. In the 
present paper, the basic ideas of the slender-body approximation are 
applied to the nonlinear transonic equation for the velocity potential 
In order to gain some insight into the essential features of slender- 
body theory at transonic speeds. The attempt has been made to present 
the results from a unified point of view which demonstrates the simi- 
larity of the slender-body solutions in the various Mach number ranges. 

The authors wish to acknowledge the valuable suggestions of 
Dr. Adolf Busemann of the Langley Laboratory during the preparation of 
this paper. 


SLENDER-BODY APPROXIMATION 


Slender-body theory deals with that class of shapes whose length 
is large compared with any lateral dimension. For such shapes at both 
subsonic and supersonic speeds, the flow in planes normal to the stream 
direction can be approximated by solutions of Laplace's equation. The 
justification is that for very slender wings or bodies the variation of 
the geometrical properties in the stream direction is small and, conse- 
quently, the rate of change of the longitudinal component of the velocity 
in the stream direction is also small. The various slender-body solu- 
tions have all been developed on the basis of the linearized potential 
equation. However, a similar development can be made on the basis of 
the nonlinear transonic equation. 

The simplest differential equation for the disturbance potential # 
which Is generally valid at transonic speeds (ref. lk, for example) is 


[1 - M 2 - (7 + DifoJjGxx + 3>yy + & zz = 0 (l) 

where x, y, and z are rectangular coordinates, M is the stream 
Mach number, and 7 is the ratio of specific heats at constant pressure 
and constant volume. With l the characteristic length and b the 
characteristic width (such as the largest lateral dimension of the con- 
figuration), the nondimens ional coordinates Xq, yq, and zq defined 

by x = Ixq, y = byq, and z = bzq and the nondimens ional potential $q 
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defined by $ = IS— <I>]_^x-j_,y-|_, are all of the order of 1 in the vicinity 
of the configuration. In this coordinate system, equation (l) becomes 


(!) 2 F - m 2 - <r + 


®1- 


•X1X1 


+ + On =0 (2) 

x yiyi 1 zizi v 


For sufficiently small values of the width parameter b /l, the terms 
involving derivatives in the stream direction can be neglected to obtain 
the result that the flow approximately satisfies Laplace’s equation 


®yy + $zz = 0 (3) 

in the crossflow plane. Equation ( 3 ) represents the slender -body approx- 
imation to equation (l) . 

The surface boundary condition is 


|i . SS(x + **) - 3a 

dn dx 


where n is the outward normal to the configuration in the crossflow 
plane. For slender configurations the surface boundary condition can 
be integrated (ref. L, for example) to give 



— dv = S ' (x) 
8n 


00 


where v is any contour enclosing the shape, S(x) is the cross- 
sectional area distribution of the shape, and the prime denotes differ- 
entiation with respect to the indicated argument. 

In the slender-body approximation, the potential satisfying equa- 
tion (l) and the surface boundary condition is represented in the neigh- 
borhood of the configuration as a solution of equation (3) plus a func- 

‘ tion of integration G(x) . Thus for r = \jy 2 + z2 < p, say, where p > b. 


$ ( x ,y,z) = 0(y,z;x) + G(x) 


(r = p) 


( 5 ) 
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where 0 is a solution of the Laplace equation in the crossflow plane 
with x appearing as a parameter introduced by the shape of the cross 
section at x. The function 0, being independent of the stream Mach 
number, can be evaluated for an incompressible flow past the shape under 
consideration. The function G(x) is determined from considerations 
involving the complete equation for transonic flow (eq. (l)) and, conse- 
quently, is dependent upon the stream Mach number and upon the shape of 
the configuration. Although the analytic expression for G(x) at tran- 
sonic speeds is not known, it will be shown that the only geometrical 
property of the configuration which influences this function is the 
cross-sectional area distribution - just as at subsonic and supersonic 
speeds . This property of G(x) is established by comparing the slender- 
body solution with the solution for the flow past a body of revolution. 

As a preliminary to these considerations it is necessary to examine the 
expression for the velocity potential in more detail. 

The flow past a slender configuration is given by the solution of 
equation (3) satisfying the boundary conditions of the problem and can 
be expressed in nondimensional terms by ~~ 

<&(x,y,z) = 




where a is the contour bounding the cross-sectional area of the con- 
figuration and/or the trailing-vortex system in the y,z plane, m is 

the unit outward normal, 0(y,z;x) = ^ S (f)' 

R = ]j(y - ti)^ + (z - and r = ^y2 + 2 2 as shown in the following 

sketch: 
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Since r is independent of the surface normal, equation ( 6 ) can be 
written as 


®(x,y,z) = ^ 

L 


'(! : ) 10S T + X(i log | da + 



elf 


(r ^ p) (7) 


where use has been made of equation (^) and where S(x) = The 


variation of C> with the azimuth angle 8 is contained entirely in the 
line integral. Two of the basic assumptions used in the derivation of 
equation (7) are that both the perturbation velocities and the perturbation- 
velocity gradients in the stream direction are small. In order to satisfy 

these assumptions B "{^) and s '”(y) must be bounded . These conditions 

imply that equation (f) applies only to shapes that are smooth and free 
from discontinuities. Moreover, an additional restriction on the asymmetry 
of the shape is sometimes required (ref. 4 ); namely, the radius of curva- 
ture of the configuration in the crossflow plane must be of the order of 
b where the shape is convex outward. 


For a body of revolution at zero incidence the contour integral in 
equation (7) vanishes and 


<&o( x >y' z ) 



(r < p 0 ) ( 8 ) 


where the subscript o is used to denote values for a body of revolution. 
Since a body of revolution is completely defined in terms of the cross- 
sectional area distribution, this is the only geometric parameter which 


enters into 


So(f)* So(f) 


is of the form 



where the 


dependence -upon the body shape is contained in s Q . Further considera- 
tion of the region of validity of the slender-body solution is necessary 


in order to show the corresponding dependence for 



Examination of equation (7) shows that the variation of the poten- 
tial with the azimuth angle becomes vanishingly small for r ^ rq, since 

the logarithm in the contour integral is of the order b/r for if « 1 . 

The magnitude of the terms neglected in equation (l) are now compared 
with those retained, in order to show that rq lies within the region 
where the slender-body solution is a valid approximation. The ratio of 
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the neglected terms 
terms for r > r^_ 




M 2 - (7 + l)M 2 $^j 
is of the order 


$xx 


to any of the remaining 




+ 0 ( 1 ) 


0 (±og 20 + 0(1) 


e 



where 0( ) denotes order of, 0(l) denotes nonsingular terms, and the 

functions g '^20 sbgL ®"(t) are cons ^ l ^ ere( ^ to te regular. From this 

ratio it can he seen that, for a given Mach number and degree of approxi- 
mation e, the region of validity of the slender-body solution (r ^ p), 
measured in terms of body widths p/b, can be made as large as desired 

by suitably restricting b/l. Consequently, for k « 1 , rj < p and 

u 

the flow field external to rj_ are nearly axisymmetric so that 


«(x,y,z) 



( r l < r < p) ( 9 ) 


In addition, for a given degree of approximation, larger values of the 
width parameter b /l are permitted at, transonic speeds than in the other 
speed ranges since the quantity 1 - M 2 is much larger at subsonic and 
supersonic speeds than at transonic speeds. _ 


In the region r > r^_, the flow about a slender configuration is 
nearly axisymmetric and $ must be identical to some $ 0 in this region. 
If <I> 0 is the potential of the associated axisymmetric flow which gives 
rise to the same velocities as <& for v > r^, then, from equations (8) 
and ( 9 ), s 0 = s and go = g. Thus, g is determined as the func-* 
tion go- Since the only geometrical property affecting g Q is s 0 , 

and since s = s 0 , the only geometrical property influencing g is s. 
Thus, the form where the dependence upon the shape 

of the configuration is contained entirely within S (^V 


In the preceding discussion the region of validity of the slender 
body approximation to was tacitly assumed to be at least as large 
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as that for <z> . This condition is certainly true since the singular 
terms in the two solutions are the same. 


A complete discussion of the validity of the slender-body approxi- 
mation at transonic speeds would require the analytic expression for 



In the absence of this information* such considerations are 


admittedly somewhat speculative. Even so* it is of interest to explore 
the nature of the approximation since some elementary considerations 
suggest that the slender-body solution will provide a reasonable approx- 
imation in regions where it might be expected to be poor - in the neigh- 
borhood of weak shock waves . Because of the nature of the slender-body 
solution* the flow is represented only in a small neighborhood of the 
configuration* and the shocks are represented as surfaces of discontinuity 
normal to the stream direction. Moreover* for slender configurations at 
transonic speeds* only near normal shock waves are to be expected. 


In the slender-body approximation the term - M 2 - (7 + 1)M‘=t&x|<D xx 

is required to be small compared with any of the other terms in the tran- 
sonic differential equation for all values of r less than p. If this 
condition is to be satisfied in the neighborhood of weak shock waves* the 
quantities 


and 







(lOa) 


( 10 b) 


must be bounded there. Since the disturbance velocities are bounded for 
shapes which satisfy the assumptions of slender-body theoiy* the quanti- 
ties in expressions ( 10 ) will be bounded at shock waves if Is 

bounded. The transonic differential equation admits of solutions having 
velocity discontinuities which are compatible with the transonic approx- 
imation to the shock-wave relations (see appendix) . Since the develop- 
ment in the appendix does not require that be singular* it seems 

reasonable to suppose that 0 xx and* hence* are "bounded in the 

vicinity of shock waves. In addition* the coefficient of &"(t) ■^ n 
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expression (10a) has a mean value of 0 for the admissible normal shock 
waves and the contour integral in expression (10b) vanishes at values 
of x/l for which the configuration is axisymmetric . 


The slender-body solutions in the various Mach number ranges are 
similar in that they are all represented by equation ( 7 ) although the 

function &(t) for the various Mach number ranges. Ward 


(ref. h) has determined the function f° r supersonic flows and 

Adams and Sears (ref. 6) have obtained a corresponding expression for 
subsonic flows. Although an analytic expression for this function at 
transonic speeds is not known., it has been established that the only 

geometric, property of the body influencing g ^0 is the area distribu- 
tion. Moreover, the transonic similarity rule for bodies of revolution 
(ref. 14 or 15 ) shows that &(t) can be expressed in the form 


g 




log 


(7 + 



where the similarity parameter is 


K = 


1 - M 2 


(7 + DM^J 


AERODYNAMIC FORCES 


Since the slender-body solutions are all represented by equation (7) 
formal expressions for the aerodynamic forces can be determined which are 
valid throughout the Mach number range. Consequently, many of the essen- 
tial features of slender-body theory at- transonic speeds can be obtained 
without resorting to detailed calculations . 


Lift 


The most significant difference between the slender-body solutions 


at subsonic, transonic, and supersonic speeds is that the function 
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differs in these various speed ranges . However, the term in the pressure 
arising from the function s(j) makes only a uniform contribution to 


the pressure at any value of x and, therefore, cannot affect the lift 
distribution or the lift. Thus, within the slender-body approximation, 
the lift distribution depends only upon the function cp and, conse- 
quently, is independent of the stream Mach number. Several investigators 
(for example, Heaslet, Lomax, and Spreiter; ref. 9) have previously noted 
that the linearized slender-body theory gave consistent results, even at 
a Mach number of 1, for planar systems . 


According to slender-body theory, the lift distribution can be 
obtained completely from solutions of Laplace's equation in the cross- 
flow plane. Since this equation is linear, the lift is proportional to 
the angle of attack even at transonic speeds . Ward has obtained an 
especially simple form for the drag due to lift in which 



where a, is the angle of attack measured from zero lift and L is the 
lift. 


Drag 

By computing the momentum change of the fluid passing through a 
cylinder enclosing the body, the drag D is determined as 


D _ £b 

1 q. 



s'(!)g'(i)d£ 



( 11 ) 


where the body extends from 5=0 to 5=1, o' denotes the contour 
of the body at the stem which in the case of wings or wing-body combina- 
tions includes the trailing-vortex sheet, q is the stream dynamic pres- 
sure, and Db is the base drag. Equation (ll) is valid throughout the 
Mach number range provided the appropriate forms of the function 

are employed. The line integral is zero for nonlifting configurations 
if the body is closed or if the body ends in a cylindrical section whose 
elements are parallel to the stream. The effect of Mach number (excluding 
the variation of base drag with Mach number) is contained in the term 

involving 6 ^ • 
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When the subsonic form of s(^J is used in equation (ll) , the 
correct result is obtained that the drag of nonlifting configurations 
is zero. By using the supersonic form of the drag varies with 

Mach number as [s 1 (l)Tj ^log ( m^ - l) . For pointed bodies, or for bodies 
which end in a cylindrical section, the supersonic slender-body theory 
indicates that the drag is independent of Mach number ■ For bodies which 
do not satisfy these conditions, the supersonic result indicates that 
the drag approaches infinity as the Mach number approaches 1. These 
results from linear theory cannot be considered satisfactory at tran- 
sonic speeds since they give a discontinuity in the drag as the Mach 
number is increased through 1 ; whereas experimental data show that the 
drag starts to increase rapidly at a subsonic Mach number and varies 
smoothly through 1. However, the few known solutions of the nonlinear 
transonic -flow equation are in good agreement with experiment in this 
regard. It would be expected, therefore, that the drag rise of slender 
shapes would be correctly approximated by equation (ll) once the tran- 
sonic form of is known. 


Transonic Area Rule 


The body shape enters into the function 



only as a function 


of the cross-sectional area distribution throughout the Mach number range. 
This property of. the slender-body solutions leads to an important result 
even though the analytic expression for gfj) is not known at transonic 


speeds. Examination of equation (ll) shows that the body cross-sectional 
shape enters into the slender-body drag expression only through the con- 
tour integral evaluated at the stern of the configuration. For a fixed 
base contour, then, the drag of nonlifting configurations depends only 
on the axial distribution of the body cross-sectional area and is inde- 
pendent of the cross-sectional shape. Thus, within the slender-body 
approximation, the drag of a nonlifting configuration is the same as 
that of the associated body of revolution having the same streamwise 
distribution of cross-sectional area provided the base contour is fixed. 
It is in this sense that an equivalent body of revolution is associated 
with a wing-body combination. This result, often referred to as the 
area rule, is especially significant at transonic speeds where larger 
values of the width parameter b /l are permitted than in other speed 
ranges . 


The property of the dependence of the drag upon the distribution of 
cross-sectional area has previously been obtained by Ward (ref. k) and 
Graham (ref. 16 ) for supersonic flow and has been observed experimentally 
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by Whitcomb (ref. 17, for example) at transonic speeds. The importance 
of this result was first noted by Whitcomb who demonstrated that the 
area rule could be used as a basis for the design of low-drag wing-body 
combinations at transonic speeds . From the preceding development, the 
transonic area rule is subject to the restrictions of slender-body theory 
with the additional condition that the base contour be fixed. 


Langley Aeronautical Laboratory, 

National Advisory Co mmi ttee for Aeronautics, 
Langley Field, 7a., January 18, 195^. 
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APPENDIX 

ON SOLUTIONS OP THE TRANSONIC DIFFERENTIAL EQUATION 
HAVING VELOCITY DISCONTINUITIES 


The transonic differential equation for the disturbance velocity 
potential (eq. (l)) can be written as 


-1 


2(7 + 1)M 1 2 dx 


[l " M2 - (7 + 1)M% X J 2 + ^ + $ zz = 0 


(Al) 


From the conservation laws, the tangential velocities across a shock 
wave are continuous and the normal velocity is discontinuous . Consider 
first the possibility that <& x is discontinuous across a surface normal 
to the x-coordinate . In order for the differential equation to admit 
such solutions, the values of <E X on each side of the discontinuity 
must give rise to the same value for the first term in equation (Al) . 
With the subscripts 1 and 2 denoting quantities immediately upstream 
and downstream, respectively, of the surface of discontinuity, this con- 
dition is satisfied by 


1 - M 2 - (7 + 1)M 2 4> 1x = _ 1 - M 2 - (7 + l)M 2 d> 2x j 


or 


1 - M 2 


(7 + 1)M 2 
2 



+ 



which is the first-order approximation to the normal-shock relations. 

Ey considering discontinuities in all three velocity components 
(i.e., oblique shock waves), the resulting expression relating the dis- 
turbance velocities on each side of the discontinuity is identical to 
the first-order approximation for the entire shock polar. Thus, the 
transonic differential equation admits of solutions having velocity dis- 
continuities which are consistent with the first-order approximation to 
the entire shock polar. Stated another way, the transonic approximation 
to the differential equation and shock relations are consistent. 



7 


MCA TN 3815 


13 


references 


1. Muni, Max M. : The Aerodynamic Forces on Airship Hulls. MCA Rep. 184, 

1924. 

2. Tsien, Hsue-Shen: Supersonic Flow Over An Inclined Body of Revolution. 

Jour. Aero. Sci., vol. 5, no. 12, Oct. 1938, pp. 480-483- 

3* Jones, Robert T. : Properties of Low-Aspect -Ratio Pointed Wings at 

Speeds Below and Above the Speed of Sound. MCA Rep. 835, 1946. 
(Supersedes MCA TN 1032.) 

4. Ward, G. N. : Supersonic Flow Past Slender Pointed Bodies. Quarterly 

Jour. Mech.. and Appl. Math., vol. II, pt. 1, Mar. 1949, pp. 75-97* 

5. Von Ka r man, Th. : The Problem of Resistance in Compressible Fluids. 

Atti del V Convegno della Fondazione Alessandro Volta, R. Accad. 
d* It alia, vol. XIV, 1936, pp. 5-59- 

6. Adams, Mac C., and Sears, W. R. : Slender -Body Theory - Review and 

Extension. Jour. Aero. Sei., vol. 20, no. 2, Feb. 1953, PP* 85-98. 

7. Lighthill, M. J. : Supersonic Flow Past Slender Bodies of Revolution 

the Slope of Whose Meridian Section Is Discontinuous. Quarterly 
Jour. Mech. and Appl. Math., vol. I, pt. 1, Mar. 1948, pp. 90-102. 

8. Keune, Friedrich: Low Aspect, Ratio Wings With Small Thickness at 

Zero Lift in Subsonic and Supersonic Flow. KTH-Aero TIN 21, Roy. 

Inst, of Tech., Div. of Aero., Stockholm, Sweden, 1952* 

9. Heaslet, Max A., Lomax, Harvard, and Spreiter, John R.: Linearized 

Compress ible -Flow Theory for Sonic Flight Speeds. MCA Rep. 956, 
1950. (Supersedes MCA TN 1824.) 

10. Spreiter, John R. : The Aerodynamic Forces on Slender Plane- and 

Cruciform-Wing and Body Combinations'. MCA Rep. 962, 1950. 
(Supersedes MCA TW's 1897 and 1662.) 

11. Lomax, Harvard, and Heaslet, Max. A.: Linearized Lifting -Surface 

Theory for Swept-Back Wings With Slender Plan Forms . MCA TN 1992, 
1949. 

Heaslet, Max. A., and Lomax, Harvard: The Calculation of Pressure 

on Slender Airplanes in Subsonic and Supersonic Flow. MCA 
Rep. 1185, 1954. (Supersedes MCA TN 2900.) 


12 . 



14 MCA TN 3815 


13* Ribner, Herbert S.: The Stability Derivatives of Low-Aspect -Ratio 

Triangular Wings at Subsonic and Supersonic Speeds . MCA TN 1423, 

19 ^ 7 - ' " 

14. Busemann, Adolf: Application of Transonic Similarity. MCA TN 2687 

1952. 

15* Oswatitsch, K. , and Berndt , S. B. : Aerodynamic Similarity at Axisym 

metric Transonic Plow Around Slender Bodies. KTH-Aero TN 15, Roy. 
Inst, of Tech., Div. of Aero., -Stockholm, Sweden, 1950. 

16. Graham, Ernest W. : Pressure and Drag on Smooth Slender Bodies in 

Linearized Flow. Rep. No. SM-13417, Douglas Aircraft Co., Inc., 
Apr. 20, 1949. " - 

17. Whitcomb, Richard T.: A Study of the Zero-Lift Drag -Rise Character- 

istics of Wing -Body Combinations Near the Speed of Sound. MCA 
RM L52H08, 1952. 


NACA - Lanpley Field, Va. 



